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Abstract In this paper, we confront the predictions of the power law cardassian model for the baryon power spectrum with the 
observations of the SDSS galaxy survey. We show that they fit only for very unusual values of the cold dark matter or baryon 
density parameters, the Hubble parameter or the spectral index of the initial power spectrum. Moreover, the best-fit Cardassian 
models turn out to be phanton models. If one wants to recover the usual values for these constants, as quoted by the WMAP 
team, the power law Cardassian model turns out to be indistinguishable from a ACDM model. 
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1. Introduction 

One of the greatest discoveries in cosmology in recent years 
is the accelerated expansion of the Universe. The first strong 
evidence that led to the conclusion that the expansion of 
the Universe appears to be accelerating, c ame in 1998 from 
observations of super novae of type la jRiess et all 1 19981 
IPerlmutter et alJ 19991). Since t hen, more recent supernovae 
obser vations jKnop et alJl2003t iTonrv et alJl2003l iRiess et al] 
120041) along with observations of the cosmic microwave 
background dSpergel et alJ 120031) and large-scale structure 



jTe gmark et all 20041) have strengthened this conclusion even 
further. 

Now that the accelerated expansion of the Universe seems 
to have been established, the task facing cosmologists is to ex- 
plain its origin. Several models have been proposed over the 
years that attempt to explain the dark energy that gives rise 
to the accelerat ed expansion. One suc h model is the power 
law Cardassian jFreese & Lewisll2002l) . This is a model that 
has no additional energy components other than baryons and 
cold dark matter. However, motivated by brane world cosmol- 
ogy JChung & FreeselfeoOOl) . the Friedmann equation is mod- 
ified by the presence of an additional energy term that is pro- 
portional to the cold dark matter density raised to a general 
power, i.e. p" t . It is this term that gives rise to the acceler- 
ated expansion or dark energy. Indeed, it can be shown that the 
Hubble function can be written in the same form as ordinary 
general relativity with a dark energy fluid defined by a con- 
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stant equation of state PdeIpde = n — 1, where pde and pde 
are respectively the pressure and de nsity of this special form 
of dark energy dAvelino et al.ll2003l) . Consequently, accelera- 
tion occurs for n < 2/3 and the power law Cardassian model 
cannot be distinguished from such dark energy models by any 
observational tests based on the Hubble function only, such as 
the redshift-luminosity distance relation inferred from super- 
novae. However, this is no longer the case for tests based on 
large- scale structure data dSandvik et alJl2004t iKoivisto et alJ 
1200.4 . since they probe more than just the background evolu- 
tion. More general Cardassian models exist where the mod- 
ification of the Friedmann equation cannot be written as a 
simple power in the matter density. Examples of such model 
are the Polytropic Cardassian Model (PC) and the Modified 
P olytropic Cardassian Mod el (MPC), which were proposed 
in iGondolo & Freesel J2003h . However, for such models it is 
more difficult to make a connection to fundamental physics, 
hence rendering the model a purely phenomenological model 
with little physical motivation. We stress that we will consider 
only the original power law Cardassian model in this work. 
Furthermore, the m odel will be considered in the so-called 
fluid interpretation ( IGondolo & Freesel 12003b . which will be 
explained in section|2] In the following, we will refer to the the 
power law Cardassian model as a phantom Cardassian when 
n < 

Constraints on Cardassian model imposed by the super- 
novae have been investigated extensively in several papers, e.g. 
dSen & Senl2003UWang et al.l2003UFrithl2004lGong & Duanl 
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120051 120061) . These studies show that the parameter space 
of the model is quite degenerate. In order to get better con- 
straints one needs to cons ider additional cosmological tests. In 
lAmarzguioui et all J2005h . the MPC model is constrained by 
comparing the matter power spectrum predicted by the model 
with that inferred from the SDSS large-scale structure data. The 
authors show that these data constrain the model severely. In 
fact, it is shown that \n\ has to be of order less than 10~ 5 , which 
renders the model virtually indistinguishable from the ACDM 
model. In the analysis leading to this conclusion the density and 
Hubble parameters w ere fixed to the first -year values quoted 
by the WMAP team JSpergel et alJl2003l) . i.e. Q m0 = 0.224, 
Q/,o = 0.046 and h = 0.72. Here, Q m o and Q,o denote the 
density parameters of dark matter and baryons respectively. In 
this paper we want to generalise this approach. We restrict our- 
selves to the power law Cardassian model, and perform a fit 
to the baryon power spectrum while allowing these cosmolog- 
ical parameters to take a range of constant values along with n. 
Furthermore, we allow also the spectral index n s of the initial 
power spectrum to take a range of possible values rather than 
being j ust unity. This make s sense in view of the last WMAP 
results JSpergel et all2006l) which give n s = 0.95 + 0.016 com- 
pared with the first-year estimate of 0.99 + 0.04. Does this 
added freedom allow the Cardassian model to deviate from the 
ACDM model? This is what we intend to answer in this paper. 



We find that all the best-fit models have negative «, which 
in the terminology defined above means that they are phantom 
models. The models that are compatible with the data fall into 
two categories: \n\ is either so small that the model is indistin- 
guishable from a ACDM model, or n can have a non-negligible 
negative value. But in the latter case, one finds that either the 
cold dark matter density, the Hubble parameter or the spectral 
index must have very unusual values for the model to differ 
substantially from a ACDM model. Thus, if one wants to re- 
cover the usual values for these constants, as quoted by the 
WMAP team, the Cardassian model turns out to be indistin- 
guishable from a ACDM model with \n\< 10 -5 . 



The structure of this paper is as follows: In section |2] we 
present the Friedmann equation of the power law Cardassian 
model and look at how first order perturbations evolve in the 
fluid interpretation of this model. In section|5]we discuss how 
to obtain the baryon power spectrum for the Cardassian model 
using the perturbed equations presented in the preceding sec- 
tion, and what initial condition to use for these. In section |4] 
we present the statistics we need to fit the predicted power 
spectrum of the model with the observed. In section |5J we 
fit the predictions of the Cardassian model with the SDSS 
jTe gmark et alJ E004I) baryonic power spectrum by keeping 
three of the parameters fi m o> ^M), h, n s fixed to their WMAP 
concordance values and allowing the remaining one to take a 
range of possible values in addition to n. 



Finally, in section[6]we summarise and conclude. 



2. Field equations 

We assume the universe to be homogeneous, isotropic and flat. 
The metric can then be written in the usual form: 

ds 2 = -dt 2 + a(t) 2 (dr 2 + r 2 (d8 2 + sin 2 0d<p 2 )) (1) 

The modified Friedmann equation for the Cardassian model is 

8ttG 



H = —^-(Pm + Pb+Pc): 



(2) 



where p m and pb are respectively the CDM and baryonic 
matter densities and p c ~ p" m . Note that we are look- 
ing at the ordinary Cardassian model and not the Modified 
Polytropic Cardassian M odel (MPC) that was considered in 
lAmarzguioui et al.N2005h . 

Energy conservation of pressureless matter implies that 
Pm,Pb ~ (1 + z) 3 and hence, p c ~ (1 + z) 3 "- This allows us 
to write the Friedmann equation as 

A 2 = (£l m0 + OmXI + z) 3 + QcoCl + z) 3 " , (3) 

no 

where we have defined the density parameters 

SnGp m o SnGpho %nGp M 



In order to derive the matter power spectrum we need to go 
one step further than Eq. Q which gives the evolution of the 
background. We need to consider first ord er perturbations. In 
order t o do this, we follow the approach of I Amarzguioui et alJ 
and consider th e fluid interpretation of the Cardassian 
model as introduced in iGondolo & Freesel J2003I) . In this ap- 
proach, the dark matter term and the term that gives rise to dark 
energy are treated as two "components" of one single fluid - the 
Cardassian fluid. We write 



P = Pm + Pc = Pm + 



_cO l-„ „ 
O PmO Pm 
"m0 



(5) 



Treating this as a perfect fluid and using energy-momentum 
conservation, we arrive at the usual continuity equation for the 
fluid p: 



p + 3H(p + p) = . 



(6) 



Substituting from Eq. © and using energy conservation of or- 
dinary matter, the effective pressure of the Cardassian fluid can 
be written as 



1 ^ d P 

P = P(Pm) = Pm~ P- 

dp,„ 



(7) 



The sound speed and the effective equation of state for the 
Cardassian fluid can now be calculated easily. We get 

2 dp d 2 p/dpf„ 
C s = = Pi 



dp dp/dp m 

«(»-l)(Q mQ + Q fc0 -l)(l+z)- 3 
«(O m0 + Om " + z)~ 3 - O,„o(l + z)" 3 " 



(8) 
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for the sound speed, and 



dp/dp,, 

I 

P 
(n 



l)(Q m0 + Q A o-l)(l+z) 



-3 



(Q,„o + Q M - 1)(1 + z)^ - Q„,o(l + z)~ 



(9) % 10000 



for the equation of state. 

In the fluid interpretation, the model we're considering con- 
tains two fluids, namely the Cardassian fluid and the baryon 
fluid. The evolution of density perturbation in these fluids is 
given by a coupled set of second order differential equations. 
The detailed derivatio n of these equations can be found in 
iLvth & Steward Jl990l) and IPadmanabhanl dl993h . Setting the 
equation of state parameter and sound speed of baryons equal 
to zero, we can write these differential equations as 

5 + H (2 - 3(2w - c 2 s )) 6 - \^H 2 (lw - 3w 2 - 6c 2 ) - (r^-fj 6 



— H (1 + w) = 

2 P+Pb 



and 

S b + 2HS b - 3 -H 2 ^l£^ 
2 p + Pb 



0. 



(10) 



(11) 



where d and 8b are the perturbations in the Cardassian and 
baryon fluids respectively. These equations are written with re- 
spect to cosmic time t . We want to solve them numerically, and 
in that respect it is more useful to write them in terms of a new 
time variable u defined by 

d 



d 

J J 

dt du 



(12) 



This implies that a = e" — 1/(1 + z) if we define the scale 
factor today as ao = 1 . Using this time parameter and mark- 
ing a derivative with respect to it with a prime, the differential 
equations take the form 



H 2 8" + (H 2 )'- + [2 - 3(2w - c 2 )] H 2 5' - Ulu - 3w 2 



6ci)H z 6 + —cid - -H 2 (l 

s ' all 1 O v 



to) 



p6 + pb&b 
P+Pb 



H ^ + {H Y±_l H 2^l£^ =0 
2 2 p+Pb 



(13) 
(14) 



3. Power spectrum 

In this work we will use the large-scale structure data of the 
SDSS team to constrain the Cardassian model. The galaxy 
power spectrum inferred from these data 1 is plotted in Fig.^ 

The power spectrum of energy component ; can be written 
as Pi(k) ~ Si(k) 2 . Thus, we can use Eqs. Jl 3i and dl4> with 
the appropriate initial conditions to predict the power spectrum 
of the baryons and/or the Cardassian flu id today. Taking the 
arguments presented in lBeca et al.ld2003l) into account, we will 
consider the power spectrum arising from the baryons only and 
not from the dark matter component of the Cardassian fluid. 

1 http://www.hep.upenn.edu/~max/sdss.html 
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Figure 1. The 22 data points of the SDSS power spectrum 

To determine the initial conditions for the coupled system 
of differential equations dl3t and H4\. we follow the a pproach 
of ISandvik et alJ (120041) and lAmarzguioui et alJ (120051) : we as- 
sume the evolution of the Cardassian and ACDM models to be 
simila r very early on, and use CMBFast dSeliak & Zaldarriagal 
1 1994 to calculate the matter transfer function at redshift z = 
200 arising from an initial spectrum with a general spectral 
index n s . This allows us to write the density perturbations at 
z = 200 as 



6(k,z = 200)~k" J2 T(z,k) 



(15) 



where T(z, k) is the transfer function obtained from CMBFast. 
Next, we use Jl 3I > and H41 to evolve the perturbations until 
today. The power spectrum of baryons today is then obtained 
as ~6 b (k,z = 0) 2 . 

4. Statistics 

To fit the 22 SDSS data points for the power spectrum we will 
use a least square test^f 2 . It is defined by 



2 (P0t,0)? bs 



X~ 



s 

1=1 



POt,0)*A0 2 



(16) 



where P(k, 0) is the power spectrum of the baryons today and 
the superscripts th and obs denote the theoretical and the ob- 
served spectra respectively. Furthermore, cr ; are the 1 <x errors 
of the SDSS data and is a normalisation factor that allows 
us to adjust the amplitude of the theoretical power spectrum. 
We want to find the values of N and the other cosmological 
parameters that minimise x 2 - To simplify this procedure, we 
can define a new quantity x 2 which minimises x 2 analytically 
with respect to the amplitude N. First, we define the following 
quantities: 



Z 



22 

A 

i=\ 

22 

!=] 

22 



(P(£,0)*) 2 



P(k, 0)f bs P(k, 0)* 



(P(k,0)° hs ) 2 

oi 



(17) 



(18) 



(19) 
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Substituting these into (I16> . we get the expression 

X 1 = AN 2 - 2BN + C (20) 

Treating this as a function of N we find that its minimum is 
reached for N = B/A. We define the value of x 2 corresponding 
to this minimising N as the new;f 2 . Thus, we write 

f = ~+C (21) 

This last expression defines a new x 2 analytically minimised 
with respect to N. This is the definition of x 2 that we will use 
to fit the data. 

5. Fit to the SDSS matter power spectrum 

In this section we perform the fitting of the theoretical 
Cardassian model to the empirical values of the SDSS mat- 
ter power spectrum. The parameters which are constrained in 
the fitting are the Cardassian parameter n, the CDM density 
parameter Q m rj> the baryon density parameter Q^o, the Hubble 
parameter h, and the spectral index n s . Ideally we would like to 
get a constraint on all of these parameters simultaneously, but 
the parameter space would be too large to handle in a simple^ 2 
test. We will therefore restrict ourselves to constraining pairs of 
these parameters with one parameter being «. 

5.1 . Constraining the n parameter only 

In this first subsect i on we mainly try to recover the results of 
lAmarzguioui et all d2005tl . In their analysis, the cosmological 
parameters were set equal to the values quoted by the WMAP 
team. However, they did not perform an explicit^ 2 fitting to the 
data. They looked at the predicted power spectrum and noted 
that it deviates strongly from the SDSS data even for extremely 
small deviation from the ACDM model. Repeating their anal- 
ysis, but doing an explicit fit to the data, we set Q m o = 0.224 
and Qm) = 0.046, implying Q. c o = 0.73. Furthermore, we set 
the value of the Hubble parameter to h = 0.72 and the spectral 
index to n s = 1. Using Eq. J2 It . we find the minimal value for 
X 2 to be 27.08 when n = -5.6 x 10~ 7 . With 22 data points and 
one free parameter, the x 2 per degree of freedom (DOF) is thus 
Xdof ~ 1-28. For n = 0, corresponding to a cosmological con- 
stant, X%of ~ 1-29. These values are not very good but they are 
comparable to the values of y 2 P ^ F that were obtaine d with the 
56 supernovae of Perlmutter ( Perlmutt er et alJl999 ) published 
in 1999, and in works aiming to constrain quintessence models 
JPi Pietro & Claeskensll2003l) . So maybe more precise data in 
the future will decrease the Xdof we 8 et m tn ^ s sect ion, just as 
was the case with the supernovae data. Figure|2]shows a plot of 
X 2 as a function of the parameter n. 

We see clearly that the SDSS data do not allow for as 
large values of n as the supernovae data: at 2 cr, n is limited 
to the range n e [-1.5 x 10~ 5 ,8.9x 10~ 6 ]. They show that 
the Cardassian model must resemble a ACDM model very 
closely to be in agreement with the large-scale structure ob- 
servations and the WMAP values for the cosmological parame- 
ters. Figure|5]compares the predictions of the Cardassian model 



with the SDSS da t a and is in agreement with the results of 
lAmarzguioui et all d2005l) . Note that negative values for n are 
not as strongly disfavoured as positive. But still only a very 
small deviation from zero in the negative direction is allowed. 
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Figure 2. x 2 as a function of n 
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Figure 3. Power spectrum predicted by the Cardassian model 
versus SDSS data for several values of n. The curves are such 
that they best fit the data. The power spectrum predicted with 
the positive value n = 0.0001 is totally excluded by the data. 

5.2. Constraining n and the CDM density parameter 

This section and the foll owing ones generalise the work of 
lAmarzguioui et all J2005h . We consider n and Q„,o to be free 
parameters and keep the other cosmological parameters fixed 
to Q/,o = 0.046, h = 0.72 and n s = 1. The minimal value for 
X 2 is then 21.86, which gives us x\>of ~ 1-09. This is obtained 
when Q„,o = 0.62 and n = -0.10. In Fig.0]we have plotted the 
X 2 along with the 1 and 2 cr confidence contours as functions 
of n and Q„,o. In Fig.[5]one can find an enlargement of this fig- 
ure for the area O,„o < 0.4. The 1 and 2<t confidence levels 
for several ranges of values of «, in particular for the small «, 
are plotted in the figure [5] At 2 cr, the allowed ranges for the 
parameters are n e [-1.1,6.3 x 10~ 6 ] and D. m0 e [0.22,0.73]. 
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Thus, we see that positive values for n are strongly disfavoured. 
But the data do allow for negative values up to order of unity 
when treating only n and Q„,o as free parameters. 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 4. x 1 an d the 1 an d 2 cr confidence contours (dashed 
lines) as functions of n and Q.,„q. We only present the values of 
X 1 smaller than 126 to clarify the colour of the plot. 




Figure 6. 1 and 2 cr confidence contours for several ranges of n. 
Only very small values of « are possible when £2„,o * 0.224, the 
value of the CDM density parameter in agreement with CMB 
data ofWMAP. 



5 

In agreement with the previous section, for the fiducial 
value £2 m o = 0.224, the Cardassian model fits the data only 
if \n\ < (9(10 -5 ), that is: it is indistinguishable from a ACDM 
model. If we want the Cardassian model to be different from 
ACDM, one has to consider some unusual values of Q m o like 
0.6. Figure[7]shows predictions of the Cardassian model for the 
power spectrum versus SDSS data for some values of Q,„o and 
n. 
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Figure 7. Some power spectra predicted by the Cardassian 
model for several values of Q„,o versus the SDSS data. 



5.3. Constraining n and the baryon density parameter 

This time, we keep Q. c q = 0.73, h = 0.72 and n$ = 1 fixed, but 
treat n and the baryon density parameter Q/,o as free parameters 
when doing the fitting. The lowest^ 2 value is now obtained for 
a negative value of Q/,o (even with a prior on Q.bo) with x 2 = 
21.55. This is clearly an unphy sical value. However degeneracy 
also allows positive values of Qm with reasonable^ 2 , e.g-X 2 = 
26.35 when fi w = 0.038 and n = -1.1 x 10~ 5 . 

Figure[8]shows a plot of the^ 2 and the 1 and 2 cr confidence 
contours as functions of f2/,o and n. An enlargement of the con- 
fidence contours is shown in figure|5]for the most relevant area 
of the parameter space. The allowed ranges for the parameters 
at 2 cr are « € [-4.3 x 10~ 5 ,5.3 x 10~ 6 ] and Q bQ e [0,0.048]. 
Once again we find that the data constrains n very strongly. Just 
as when we constrained n alone, we find that the allowed range 
is smaller for positive than negative values. But the constraint 
is still so strong that n has to satisfy \n\ < 10~ 5 in order for the 
Cardassian model to be in agreement with the SDSS data. 

Figure [K)] shows predictions of the Cardassian model for 
the power spectrum for some choices for n and Q^o versus the 
SDSS data. We see from Fig. |9] that the data seem to prefer 
values for Sl^o that are smaller than the WMAP value. But still, 
the allowed values of n are so small that it leaves the Cardassian 
model virtually indistinguishable from a ACDM model. 
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Figure 8. x 1 an d the 1 and 2cr confidence contours (dashed 
lines) as functions of Q/,o and n. Only the small value of « allow 
to fit the SDSS data. 
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Figure 9. The 1 and 2<x confidence contours in the plane 
(n,Q M ,). 

5.4. Constraining n and the Hubble parameter h 

Here, we treat n and h as free parameters in the fitting, and 
keep Q,„o = 0.224, Q/,o = 0.046 and n s = 1 fixed. The best 
fit is now x 1 = 21.02 for h = 1.11 and n = -9.8 x 10~ 5 with 
Xdof ~ 1-05. Plots of the x 1 an d the 1 an d 2cr confidence 
contours as functions of h and n can be found in figure ^2 A 
couple of enlargements of the confidence contours are plotted 
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Figure 10. Some power spectra predicted by the Cardassian 
model for several values of Q&o versus the SDSS data. 



in figure El The parameters ranges allowed at 2cr are n < 
1.8 x 10 5 and h > 0.76. We have not calculated the lower 
and upper limits for respectively n and h since they are beyond 
h = 1 .40 (and then below n — -8 x 10~ 4 ). Thus, the data appear 
to place a tight limit on positive values for «, but not on negative 
values, at least if one doesn't put a prior on h. Clearly, as the 




Figure 11. x 1 an d the 1 and 2 sigma confidence contours (dot- 
ted line) as functions of Ho and n. 



confidence contours show us, the SDSS data are not able to 
constrain the Hubble parameter very well. 

Note that the HST Key Project jFreedman etail l200lh 
quotes a value for the Hubble parameter of h = 0.72 + 0.8, 
whereas the WMAP team quote a value of h = 0.72 + 0.05 
when using WMAP data alone. If we demand h to be around 
0.7 in order to be in agreement with these measurements, we 
find that \n\ must be of order less than <9(10~ 5 ). Again, this im- 
plies that the Cardassian model is indistinguishable from the 
ACDM model. The only way to make it substantially different 
is to accept values for h larger than 1 .0, which is quite unreal- 
istic with respect to the measurements quoted above. 

Fig. 1131 show plots of the power spectra predicted by the 
Cardassian model for some choices for the parameters n and h. 
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50000- n=-0. 00009, H = 110.93, x 2 =21.08 
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Figure 13. Some power spectra predicted by the Cardassian 
model for several values of Ho versus the SDSS data. 

5.5. Constraining n and the spectral index n s 

In all the previous sections we have chosen the Harrison- 
Zeldovich power spectrum P - k" s with n s — 1 as the initial 
power spectrum. In this section we will consider n s as a free 
parameter along with «. The other cosmological parameters are 
kept fixed with Q w = 0.046, Q m0 = 0.224 and h = 0.72. The 
best fit is then^ 2 = 21.25 with Xdof ~ when n s = 1.23 
and n = -4.0 x 10~ 5 . The^ 2 and the 1 and 2 <x confidence con- 
tours are plotted in Fig.^Jas functions of n s and n. At 2cr the 
allowed ranges for the parameters are n e [-0.3,4.2 x 10" 6 ] 
and n s e [1, 1.73]. Just as we saw in the constraints in sec- 
tions |^] and positive values for n are tightly constrained, 
whereas relatively large negative values are allowed. 

Thus, it is evident that the SDSS data do not tightly con- 
strain n s . Indeed, we find that for n s 1.6, large values of 
n < are in agreement with the SDSS data. Thus, if we ac- 
cept such large values for the spectral index, it is then possible 
to construct a Cardassian model that differs considerably from 
ACDM and still agrees with the SDSS data. However, such 
values of n s are not compatible with the value measured by 
WMAP. For the special value n s — 1, i.e. a Harrison-Zeldovich 
spectrum, we find that \n\ has to be of order less than 0(1 -5 ) 
in order to be in agreement with the data. 

Once again, if we do not modify the usual Harrison- 
Zeldovich power spectrum with n s = 1, the Cardassian model 
is equivalent to a ACDM model. Fig. ^] shows plots of the 
predicted power spectrum for the Cardassian model for several 
values of n and n s along with the SDSS data. 

6. Conclusions 

The purpose of this work has been to constrain the parame- 
ters of the power law Cardassian model in the fluid interpreta- 
tion by using large scale structure data from the SPSS galaxy 
survey. We generalised the work of I Amarzguioui et alJ J2005t) 
by treating additional cosmological parameters as free in ad- 
dition to n itself under a fitting to the galaxy power spectrum. 
In our analysis, we looked at five different fits to the SDSS 
galaxy power spectrum. In the first fit, repeating the analysis of 
Amarzguioui et al. d20Q5l> . we fixed the four cosmological pa- 




0.02 0.05 0.1 0.2 

k (h Mpc" 1 ) 



Figure 15. Some power spectra predicted by the Cardassian 
model for several values of n s and n versus the SDSS data. 



rameters Q mU (CDM density), Q&o (baryon density), h (Hubble 
parameter) and n s (spectral index of the initial spectrum), and 
treated only the Cardassian parameter n as a free parameter. 
One finds then that \n\ has to be of order less that 10~ 5 in or- 
der to be in ac cordance with t he da ta. A similar conclusion 
was reached by iKoivisto et alJ d2005l) for the MPC model us- 
ing the CMB power sp e ctrum . These results are also similar to 
those of ISandvik et alJ J2004I) . who showed that a generalised 
Chaplygin gas has to resemble the ACDM model closely to 
avoid oscillations or blow ups of the matter power spectrum. 

In the remaining four fits, we allowed one additional pa- 
rameter to take values freely along with «, while keeping the 
remaining ones fixed to their WMAP best fit value. First, con- 
straining f2,„o simultaneously with n, we find that positive val- 
ues for \n\ are tightly constrained, but not the negative. While 
the 2 <x upper limit for n is of order 10~ 6 , the lower limit al- 
lows for negative values of order unity. But for such negative 
values the Universe would have to contain an unsual amount of 
cold dark matter. If we demand that the dark matter content be 
within the bounds given by WMAP, we again find that \n\ has 
to be so small that the model becomes indistinguishable from 
ACDM. 

Next, varying Q/,o and n freely, we again find that latter is 
strongly constrained both in the negative and positive direc- 
tions. More specifically, \n\ has to be of order smaller than 10~ 5 
to be in agreement with the data. Next, varying h and n freely, 
we again find that positive values for n are tightly constrained, 
but negative values are much less constrained. However, it turns 
out h has to be larger than ~ 1 .4 in order for n to be smaller 
than -10" 3 . This is an unrealistically high value for the Hubble 
parameter. Demanding a more realistic value for h will again 
render the Cardassian model indistinguishable from ACDM. 

Finally, allowing the spectral index n s and n to vary freely, 
we find the same behaviour as in the last fit. Positive values 
for n are tightly constrained, while negative values are not. But 
adding a reasonable prior on n s will restrict the allowed nega- 
tive values severely. In fact, n s has to be larger than ~ 1 .5 in 
order for n to be of order less than -10 -1 and still agree with 
the SDSS data. A summary of our results showing the best fit 
models can be found in tabled 
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Hence, unless one is ready to accept very unusual values for 
the cosmological parameters considered here, the SDSS data 
force the power law Cardassian model to be virtually indistin- 
guishable from an ordinary ACDM model. It thus seems dif- 
ficult to unify dark energy and dark matter via the power law 
Cardassian model. 

Does this spell the end for the power law Cardassian model 
- at least in the fluid interpretation? A more general treatment 
is probably needed to conclude this absolutely. For example, 
one would ideally like to write a code like CMBFast for the 
Cardassian model instead of using the ACDM model to gener- 
ate the initial perturbations. However, this work confirms that 
the SDSS data look like a promising way to rule out this model. 
It shows that the galaxy power spectrum impose much tighter 
constraints on the Cardassian model than the supernova data. 

Finally, we wish to stress that we have only considered the 
Cardassian model in the fluid interpretation. If an alternative 
interpretation were to be found, where perturbative calculation 
could be performed, a similar analysis might very well produce 
a different result. 



Free parameters Priors 



XboF Best fit 



n n s = l, n ra0 = 0.224, 1.28 n = -5.6 x 1(P 

Qim = 0.046, h = 0.72 

n, C2 m0 n s = 1, Q.,,a = 0.046, 1.09 C2 m0 = 0.62, 

/? = 0.72 » = -0.1 

n, n w n s = 1, Q c0 = 0.73, 1.08 n w = -0.0022, 

h = 0.72 n = -2.0 x 10~ 5 

n,h n s = 1, n m0 = 0.224, 1.05 h = 1.11, 

Qjo = 0.046 n = -9.0 x 10~ 5 
n,n s £2,,,o = 0.224, 1.06 n„ = 1.23, 
= 0-046, h = 0.72 n = -4.0 x \0 S 

Table 1. Best values of the cardassian parameters fitting the 
SDSS baryon spectrum 
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Figure 12. x 1 as a function of Ho and n. 




Figure 14. x 2 and the 1 and 2 sigma confidence contours (dotted lines) as functions of n s and n. Note that in the second plot, the 
lcr confidence contour is not visible at this scale. 



